Linear AC Response of Diffusive SNS Junctions 
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We explore the behavior of the ac admittance of superconductor-normal metal-superconductor 
(SNS) junctions as the phase difference of the order parameters between the superconductors is 
varied. We find three characteristic regimes, defined by comparing the driving frequency uj to 
the inelastic scattering rate F and the Thouless energy Et of the junction (typically hV <^ Et)- 
Only in the first regime uj <^ V the usual picture of the kinetic inductance holds. We show that 
the ac admittance can be used to directly access some of the characteristic quantities of the SNS 
junctions, in particular the phase dependent energy minigap and the typically phase dependent 
inelastic scattering rate. Our results partially explain the recent measurements of the linear response 
properties of SNS Superconducting Quantum Interference Devices (SQUIDs) and predict a number 
of new effects. 



The frequency dependent susceptibility typically re- 
veals information about the internal dynamics of the 
studied systems. In the electronic case, this susceptibility 
is more often measured as admittance, whose frequency 
dependence in semiclassical models for bulk metals is due 
to scattering and appears for frequencies exceeding some 
tens of THz [1 , and in wires it is dictated by stray ca- 
pacitances and geometric inductance. Understanding the 
frequency-dependent response is moreover of importance 
for high-frequency devices. 

In superconductors or superconducting tunnel junc- 
tions (SIS), the admittance at low frequencies is domi- 
nated by the superconducting kinetic inductance [2] . For 
SIS junctions, this Josephson inductance is related to the 
supercurrent Is{^) through the junction, which depends 
on the superconducting phase difference (p: 

Op 

LAv>r' = -jrd^Isi^) . (1) 

This relation allows characterizing the current-phase re- 
lation of Josephson junctions via measurements of their 
ac admittance [3 . The remaining dissipative part of the 
admittance is due to quasiparticles, is proportional to 
exp[— |A — huj\/{kBT)] [4 , and is important only for the 
high frequencies of the order of the superconducting gap 
A or temperatures T close to the critical temperature. 

For other types of Josephson junctions than SIS the 
admittance, however, can deviate from the above simple 
picture. In this Letter we show how combining normal 
metal wires (N) and superconductors (S) into SNS weak 
links results into an admittance that entails characteris- 
tics of the inelastic scattering rates and the inverse dif- 
fusion times through the structure. At frequencies of the 



order or larger than the inelastic scattering rate F, the 
simple kinetic inductance picture has to be revised to in- 
clude non-adiabatic effects associated with the dynamics 
of the electron distribution. The dissipative response, de- 
scribing microwave absorption, is moreover finite for tem- 
peratures or frequencies exceeding the phase-dependent 
minigap in the spectrum of excitations inside the junc- 
tion. It probes the density of states in the junction, and 
is related to the physics of stimulation and suppression 
of the supercurrent O [6] . 

Here, we study diffusive SNS junctions, whose length 
L is longer than the superconducting coherence length 
^0 = A, where D is the diffusion constant. The 

proximity effect induces a gap in the density of states 
inside the normal metal, Eg{(p = 0) = 3.12Et <C A, 
where Et = hD/L^ is the inverse diffusion time. The 
minigap depends on the phase difference approximately 
as cos((/p/2) and vanishes for cp = tt. The coupling to the 
electric field is modeled by assuming an ac bias voltage, 
V{t) = 5Vsm{ujt)^ which induces an oscillating phase 
difference (p{t) = (p -\- 6(j)cos{ujt) across the junction. 

To find quantitative results, we describe the SNS junc- 
tion dynamics with the Keldysh-Usadel equation [7] , used 
also in Ref. [6 . In this approach, physical quantities 
are obtained from the Retarded, Advanced and Keldysh 
Green's functions g^^^^^ {E^ E')^ which depend on two 
energy arguments. These functions are matrices in the 
Nambu (electron-hole) space, and the Keldysh (K) part 
can be parameterized in terms of an electron distribu- 
tion function mditiix h{E,E') = hL{E,E')^hT{E,E')%\ 
_ -Rj^_j^-A ^ where matrix products involve also con- 
volutions, {BC){E,E') = j^^^B{E,E^)C{E^,E'). 
In the presence of the harmonic drive, these functions 



can be written in a matrix representation gn^rn{E)=g{E-\- 
nhuj^E + mhu) [8 that reduces convolutions to matrix 
products. 

The ac admittance is Y{uj) = j^j^^I {u) ^ for a hnear- 
response drive 5(j) <C 27r. The ac current harmonic is 
= ^ dETirsJoiiE), where aN and S are the 
normal-state conductivity and the cross section of the 
junction, and the Keldysh current is = g^dxQ^ + 
g^OxQ^' Here x is the position in the junction, and 
dxB = dxB — i[Afs^ B] the gauge-covariant gradient con- 
taining the vector potential A{E, E') = A{E)5{E - E'). 

In the above approach, the admittance splits naturally 
into three gauge-invariant parts, F = Fgc + ^dy + ^qp, 
where (hereafter, h = e = ks = I) 
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^ J JETr{Mj^,hu-hooj^,)}, (2a) 
1^ dE Tr{fsCAi - hodAn , (2b) 

|^£d^TV{(l-f3^f,f3^„-0) (2c) 
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Here j^/^ = g^^^dxQ^^^ describe spectral (su- 
per)currents, hoo{E) = hii{E — uj) = tanh(£^/2T) 
the equilibrium electron distribution, and /iqi its time- 
dependent part. The contributions describe (a) the ac 
supercurrent, (b) effect of the dynamic variation of the 
populations of the Andreev levels, and (c) the quasiparti- 
cle current driven directly by the field. Below, we mostly 
work in a gauge in which the electric field is contained in 
the vector potential, A{t) = Aocos{ujt)^ Aq = —5<p/{2L). 
Requiring charge neutrality leads to a finite position de- 
pendent scalar potential, but our numerics indicate that 
this can be disregarded for u < IOEt- 

In the following, our aim is to relate the contributions 
([2| to quantities that can be calculated in the absence 
of the ac drive [9]. The regimes where the different con- 
tributions are relevant depend on the particular values 
of the phase difference, frequency and temperature. The 
main results are summarized in Fig. [ij which shows the 
phase dependence of the admittance in different regimes 
of frequencies and temperatures. 

Low frequency: For frequencies satisfying F <C <C 
Eg{ip)^ the superconducting correlations follow the time- 
dependent phase difference, but the electron distribu- 
tion is driven out of equilibrium. In this limit, the 
time-dependent supercurrent (2a) yields Eq. Q, i.e., 
Fgc :^ —2idipIs{^)/oj. Since the supercurrent Is decays 
exponentially as the temperature increases, this contri- 
bution to reactance becomes unimportant at high tem- 
peratures (T > 10£^t), unless the frequency is very low. 

The second major contribution to reactance comes 
from a dynamic variation of the population, as given 
by Idy For this, the time dependent component /loi = 
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FIG. 1. (Color online): Phase dependence of the admittance 
Y in units of the normal state conductance Gn (dissipative 
part: solid, reactive part: dashed), in different regimes of 
interest. The "low" temperature results have been calculated 
for T — Et, the "high" temperature results for T — IQEt- 
The adiabatic frequency is cj = ^t/200 < F = Et/2^, the 
low UJ = Et/^ ^ F, and the high frequency is cj = IO^t- 
Dotted lines show the contribution from dominant parts, 
[Eqs. ([l]), 11}], ^dy [Eq. ([s])], Fqp [Eqs. (|9])], indicated in 
the lower left corners. The adiabatic frequency results are 
obtained from analytical approximations, the others from the 
full numerics. 



h']^ -\-rsh'rp of the distribution function needs to be solved 
from a kinetic equation. Assuming again simple time de- 
pendence for the superconducting correlations, the first 
harmonic of the Usadel kinetic equation reads (cf. 0) 



Dd, ■ iVLd^h'^ - Td^h'j, + jsK) = (3a) 
\js + 5,r)(/iii - /loo) + i(w - 2iV)Nh!^ , 
Dd^ ■ {Vrd^h'T + m/il + jsh'L) = 0. (3b) 



Here Pl/t are the spectral heat /charge diffusion coef- 
ficients, T is an anomalous kinetic coefficient and N 
is the local density of states. These quantities are 
related to the equilibrium Retarded Green's function, 
g^{E) = g{E)Ts + f{E)T^ - f{E)T^, e.g., js = Im j^;, 
JE = h TrfaJol), N = Re^, and Vt = ^ + + 72+ 
l/T/2]. 

At low frequencies, Eqs. (|3| should be solved with 
Andreev reflection boundary conditions amounting to 



and dr • h'r 



at the two NS interfaces. The re- 



sulting function hip is in the vector potential gauge flnite 
but small, and as a flrst approximation we can disregard 
it. Moreover, gradients of h'j^ are small due to Andreev 
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FIG. 2. (Color online): (a) Function Q{(p,T) (see Eq. ^) 
describing the phase-dependent dynamic contribution to the 
reactive response and (b) function P{ip,T) (Eq. ([7|) de- 
scribing the phase dependence of the dissipative part of 
the admittance at low frequencies. From top to bottom 
in (a) and marked with symbols in (b): T/Et — 16(A 
),8(v),4(n),2(o), l(+),0.5(x). The dotted lines represent 
the analytic high-temperature approximations to which Q and 
P tend for T ^ Et- The inset of (b) shows the temperature 
dependence of Go(T). 



8.9(1 saw((/?)| — 0.32 sin^((^/2)) I sin((/?)| where saw((^) = 
{(f -\- TT mod 27r) — tt. At low temperatures, Q{(p^T) is 
suppressed for phases at which T < Eg{(p). Note that 
Q is positive definite, has almost a double periodicity 
compared to the kinetic inductance term, and has a min- 
imum around (/p = 0, where the kinetic inductance is at 
maximum. 

The dissipative part of the impedance originates from 
two additional sources: the quasiparticle part Yqp, and, 
importantly at low temperatures, a part of the AC super- 
current oscillating in phase with the voltage. The quasi- 
particle contribution is easiest to derive in a gauge where 
the vector potential vanishes. Then, hip in Eqs. ([3| has 
the boundary conditions h'rp{x = ±i^/2) = ±^(/iii — 
/loo), and solving Eq. (|3| assuming h'j^ ^ yields 
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Go(T) 



4rcosh2[^/(2T)] 



(6) 
(7) 



reflection, and we get a fairly good estimate for the av- 
erage h^r by averaging Eqs. (|3j) over the normal-metal 
junction, defining (•) = dx ■ /L. As a result, we get 



2 {uj-2ir){N)' 



Substituting this to Eq. (2b) and assuming uj < Et, 



dy 



-iGn Ej 
-2iTY 

- 2iT 



(4) 



(5) 
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J — C 



dE- 



AT{N) cosh^[£'/(2T)] 



This is similar to a correction to the dc conductance 
described by Lempitskii [10]. Its origin can be under- 
stood as follows [llj: the current is carried by a dense 
spectrum {en{^)} of discrete bound states with popula- 



- ^^ni^^^n) fni^n) • With ac bias one finds 
\enfn/i^^Y.n(d^enfdJn/i^, whcrc the 



tions /n, j(^) 

first term is equivalent to Fgc and the second one to Y^ 



The dynamic contribution is purely dissipative and 
constant for <C T, contains both reactive and dissipa- 
tive components for ^ F, and becomes purely reactive 
and decays for F, as visible in Fig. [T] In general, 

F also depends on the phase difference [12], so that the 
phase-dependent response at frequencies of the order of 
the inelastic scattering rate may be quite complicated. 
On the other hand, Eq. ([5| offers a way to probe such 
phase-dependent scattering rates via an admittance mea- 
surement. 

The function Q{(f,T) is shown in Fig. |2|a) for 
different temperatures. At T ^ Et the re- 
sponse can be fitted with the function Q{(p) ~ 



where we assumed uj < Et, and defined Go{T) as the 
(p = value [P(0,T) = 0]. The kernel is Kq ^ 
(D-i)-i - Red^jE ^ + + [13. It de- 

scribes the spectrum of excitations in the junction avail- 
able for receiving energy — this has a minigap Eg{(p)^ so 
that for uj^T < Eg{(p) all dissipation vanishes. Note that 
the appearance of the minigap is related to the presence 
of the AC supercurrent contribution: Fqp > Gat and is 
similar to the usual proximity-enhanced conductance. 

The functions P{(p^T) and Gq{T) are shown in 
Fig. [2|b). At low temperatures, the temperature and 
phase dependence shows a clear signature of the pres- 
ence of a minigap in the density of states; note that this 
also applies to the Idy contribution. The dissipation is 
concentrated at phase differences close to tt where the 
minigap is small. In the high-temperature limit, the dissi- 
pative term consists of a phase independent contribution 
Go{T) and the phase-dependent part has a simple phase 
and temperature dependence, — 0.23 saw((/p)^GAr£^T/P- 

High frequency: When the frequency becomes of the or- 
der of the Thouless energy Et, the above semi-adiabatic 
expressions break down as the spectral quantities become 
frequency dependent. We can however construct approx- 
imations also in this limit (for full expressions, see Ap- 
pendix). The supercurrent contribution is fairly approx- 
imated by 



Yac 



4uj 
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d^[jE{E)+jE{E + uj)]*hoo{E) 



■d^[jE{E)+jE{E + uj)]hi^{E) 



(8) 
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FIG. 3. (Color online): Amplitude 5Y 

maXc^[Re / Im]y((/?) — minc^[Re / Im]y((/?) of the phase- 
dependence in the nonadiabatic dissipative (Re) and reactive 
(Im) admittance of a long SNS junction, obtained by solving 
the time-dependent Usadel equations numerically. The lines 
correspond to temperatures T / Et = 0.5,1,2,4,8,16 (top to 
bottom) and were calculated with V — 0.05^t- 



and the quasiparticle part of the impedance by 



/oo 
dEK{E,oj) 
-OO 



hn{E)-hooiE) 
2uj 



(9) 



K{E,uj) = -{fi.+g{E)g{E + ujy (10) 
+ \f{E)f{E + oj)* + \f{E)~f{E + uY] - 

The dynamic contribution l^y can be neglected for uo ^ 
The above is compared to full numerical solutions 
in the top panels in Fig. [T] 

At high temperatures, Fgc is exponentially suppressed 
even for high frequencies, similarly as the equilibrium su- 
percurrent. Consequently, Yqp dominates for cj, T ^ Et 
(see Fig. [T]). At low temperatures, both the supercurrent 
and quasiparticle contributions are important. 

The amplitude of the phase dependence in the ad- 
mittance is illustrated in Fig. |3) up to high frequencies. 
As the frequency increases, the proximity-induced phase 
dependence in the reactive and the dissipative compo- 
nents decays, as the admittance approaches the constant 
normal-state value, Y{uj^(p) Gn- Similarly, the phase 
dependence vanishes as the temperature increases. At 
low frequencies, on the other hand, the reactive compo- 
nent diverges due to the Josephson inductance, and the 
dissipative component is dominated by Idy 

In a recent experiment probing directly the ac admit- 
tance of SNS junctions [11^, it was found that at frequen- 
cies r <C cj <C Et (semi-adiabatic limit), the reactive re- 
sponse follows closely our prediction consisting of the sum 
of the Josephson inductance and the dynamical correc- 
tion ([5|. However, both the dissipative contribution and 
the dependence at high frequencies u > Et are different: 
in [11 , the dissipative contribution is directly related to 
the reactive contribution, and moreover the amplitude of 
phase oscillations in susceptibility Sx = iojSY decays as 
frequency increases. The characteristic frequency scale 



for this was found to be temperature independent and of 
the order of the Thouless energy. In contrast, in this fre- 
quency range our Usadel model predicts Sx ^ const. — 
however, it is for example possible that the simple relax- 
ation time approximation does not include all interaction 
mechanisms playing a role in the experiment. 

Finally, we remark that our work amounts essentially 
to deriving the parameters for the Resistively Shunted 
Junction (RSJ) model of SNS junctions [2 : there, the 
Josephson inductance resulting from the supercurrent 
term should be modified to include the nonadiabatic cor- 
rection, Eq. ([5|, and the shunt resistor describing the 
dissipation in the junction should be replaced by the dis- 
sipative terms presented in Eqs. (5j7). 

In conclusion, we have described the frequency- 
dependent admittance of diffusive superconductor- 
normal metal-superconductor junctions and shown how 
the simple adiabatic Josephson inductance picture is 
modified once the frequency is increased. Besides study- 
ing the dynamics of the system, the detailed frequency 
dependence can be used to study directly the inelastic 
scattering rates. Our results are also relevant for devices 
utilizing high-frequency properties of SNS junctions, such 
as those used in metrology and radiation detection. 
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admittance can be obtained with quantitative accuracy if 
the solution at equilibrium is known. The approximation 
procedure is outlined in the main text. Below, we show 
intermediate results before taking the uj ^ limit, and 
demonstrate that the results agree with the full numerical 
approach which solves the ac Usadel equation exactly (see 
the Supplementary information of Ref. 6 for details of the 
numerics). 



On approximations 

Although one cannot solve the time-dependent linear- Following the procedure outlined in the main text, the 

response Usadel equations analytically in closed form, the following approximations can be derived: 

I 



/oo nOO 

dE Y{E) = aS dE + fdy(-E) + Yqp{E)] , (11a) 
-OO J — OO 

^^^(^) ^ zi (d4ME)±Ml±^,^E + - Mm±ME±^niE)) , (lib) 

2uj V 2 2 / 

- -^^#f^^^^4^(M^ + -) - W), (lie) 

%v{E) ^ ^([1 + 9{E)g{E + uj)* + \f{E)f{E + w)* + \hE)f{E + w)*]-i)-'(/i(i? + w) - h{E)) . (lid) 



Above, h{E) = tanh[£^/(2T)], and the energy-dependent 
quantities jE{E) = ^ Tr[f3^^V^^] and g^{E) = 



g{E)h^f{E) 



ri+if2 
2 



f{E) ^^ J'^'^ are obtained from the 



equilibrium Usadel equations f7|, which are much faster 
to solve than the full ac equations. 

Figures[4j[5) and [6] illustrate that the above approxima- 
tions reproduce all qualitative features visible in the fully 
numerical solution, and are quantitatively accurate in a 
large part of the parameter regime we are interested in. 
In particular, the deviations between the two approaches 



are miniscule for frequencies lower than Et- At large fre- 
quencies there are clear quantitative differences, but the 
qualitative phase and frequency dependence is the same, 
and the resulting admittances are of the same order of 
magnitude. This demonstrates that to a fair accuracy the 



ac admittance can be well described by using Eqs. (11) 



and standard solvers for the equilibrium Usadel equation. 



Besides that, Eq. (11) provides a possibility for making 



analytical estimates for the admittance contributions. 
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FIG. 4. Left: Fig. 1 in the main text. Right: Fig. 1 with all data computed from Eqs. (11). 




FIG. 6. As Fig. 3 in the main tex t, b ut expressed as a color plot. Left: numerically computed admittance oscillations. Right: 
admittance oscillations from Eqs. (11). 



